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The calculation of the local density of states (LDOS) in lossy materials has long been disputed due to the di- 
vergence of the homogeneous Green function with equal space arguments. For arbitrary shaped lossy structures, 
such as those of interest in nanoplasmonics, this problem is particular challenging. A non-divergent LDOS ob- 
tained in numerical methods like the finite-difference time-domain (FDTD), at first sight, appears to be wrong. 
Here we show that FDTD is not only an ideal choice for obtaining the regularized LDOS, but it can address the 
local field problem for any lossy inhomogeneous material. We exemplify the case of a finite-size photon emitter 
embedded within and outside a lossy metal nanoparticle, and show excellent agreement with analytical results. 
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The spontaneous emission rate is proportional to the 
imaginary part of the photon Green function (GF) with 
equal space arguments, i.e., F oc Im[G (r,r; w)], and can 
be decomposed into a homogeneous contribution and a 
scattered contribution, G = G*^"™ -I- G^'^'***. In the case 
of a lossless homogeneous material, then 



Im[Gr^(r,r;w)]=fc„V/67r, 



(1) 



while the real part diverges (here fco is the vacuum 
wavevector, n is the refractive index, and fi is the per- 
meability of the homogeneous space). As soon as loss 
is introduced (e.g., Im[n],Im[/z] ^ 0), then the real and 
imaginary parts of the GF are mixed, causing both to di- 
verge. This problem is well known and has been discussed 
in the context of quantum optics for decades [1,2]. 

One commonly proposed solution to this LDOS diver- 
gence has been to consider the local environment around 
the (photon) emitter to be a lossles.s cavity [1,3], thus 
modelling the effects of the local field as seen by the atom 
instead of the macroscopic field. This circumvents the is- 
sues arising due to the divergence of G(r,r), while still 
including the contribution of the outer material response 
via G*"^''**. For larger photon emitters such as quantum 
dots (QDs) and macromolecules, the dipole model in a 
fictitious cavity may not be the best approach. Moreover, 
most of these local-field models are restricted to "spheri- 
cal cavities," and are usually (but not always) applied to 
homogeneous structures. For more general nanophotonic 
structures, a common numerical method of choice is the 
finite-difference time-domain (FDTD) technique [4]. In a 
lossy structure, FDTD obtains a finite G(r,r) for both 
the real and imaginary parts, and this raises the ques- 
tion about possible numerical problems near and within 
lossy structures. This question of FDTD applicability in 
lossy nanostructures becomes even more pertinent as re- 
searchers begin to build hybrid plasmonic/photonic sys- 
tems [5], which contain QDs and metal particles. Many 
colloidal QDs are also known to be described in terms of 



a finite-size dipole within a complex index [6]. 

In this Letter, we demonstrate that FDTD can be ef- 
ficiently applied to compute the regularized GF in any 
inhomogeneous lossy medium. We compare with known 
exact solutions and show excellent agreement. We also 
address the local-field problem, and again show that 
FDTD can calculate the local-field G for any arbitrary 
structure. As a representative example, we consider a 
spherical metal nanoparticle (MNP) in air with a radius 
a = 20 nm where the permittivity of the MNP is given 
via the Drude model, e = Sr — lo^jiyP' + 47^), with pa- 
rameters similar to silver: e^ — 6, LJp/2'K — 7.89 eV, and 
7/27r = 51 meV. We define a projected local density of 
states (LDOS) along the ith axis as. 



/Oj(r;w) = 



lu){Gii{Y,v;uj)] 
Im[Gyf'=(r,r;w)]' 



(2) 



where we have normalized by the vacuum GF, so that 
this function directly gives the Purcell factor [7]. 

In Fig. 1 we show schematics of the geometries we 
will investigate. Figure 1(a) depicts a homogeneous space 
where we have divided the region onto a grid. For this 
scenario, we calculate the GF via two methods: {i) using 
FDTD and varying the FDTD grid size (Yee cell [4]), 
and {ii) analytically integrating the homogeneous GF 
over a cubic integration volume [8] . Both of these meth- 
ods act to regularize the GF over a finite size. The need 
to obtain regularized GFs is known to produce physically 
meaningful results, even in the context of point scatters 
in free space [9]. Figure 1 (b) shows a similar calculation 
geometry as in (a) except we now consider a MNP, where 
the total GF can be broken into homogeneous parts [as 
in Fig. 1(a)] plus scattered parts [10] — and compared 
with the total GF as calculated using FDTD. Figure 1(c) 
shows the real cavity model where a small cavity of loss- 
less dielectric (vacuum in this example) is placed inside 
a homogeneous space [1] which can be decomposed into 
homogeneous parts, Eq. (1), and scattering parts using 
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Fig. 1: Schematics of the regularized GF integrated over 
a finite volume, for (a) a homogeneous material, and (b) 
a MNP. Schematics of the real cavity GF are shown in- 
side (c) a homogeneous material, and (d) a MNP. (e) 
Py{r) as a function of height above and inside a MNP. 
FDTD calculations (log scale) are show for 1-nm (blue 
circles) and 2-nm (orange crosses) grids. Outside (above) 
the MNP, these are compared with the analytic GF (solid 
green line). The solid vertical black line shows the edge 
of the MNP, and the dashed vertical black lines show the 
cases considered in Fig. 2. 

analytical methods [10]; or the total GF can be calcu- 
lated using FDTD. Figure 1 (d) shows the calculation 
geometry of the real cavity model inside a MNP of ra- 
dius a, where the emitter cell has a different (and real) 
refractive index; this latter geometry can be decomposed 
into a spherical multilayer problem [3] , or solved directly 
via FDTD. We stress that we can apply the FDTD to 
any arbitrary shaped structure, but we choose a spherical 
geometry to compare with known analytical solutions. 

Figure 1(e) shows the peak LDOS along the y direc- 
tion, /9j^(r, cjpoak), for the MNP as a function of height 
along the z direction. We plot FDTD results using blue 
circles (orange crosses), corresponding to a gridding of 
a/20 (a/10), both of which corresponds to the scenario 
shown in Fig. 1(b). We use FDTD Solutions [11] and 
employ conformal meshing [12] which is used to reduce 
staircasing effects due to the rectangular grid. We see 
that inside the MNP (z/a < 1) the LDOS peak is al- 
most constant up until the very edge of the MNP, in- 
dicating that the homogeneous contribution dominates 
this region. Additionally, the two grid sizes give differ- 
ent values for the LDOS only inside the MNP, despite 
ensuring we have more than 20 FDTD grid points-per- 
wavelength. Different results for reduced grids typically 
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Fig. 2; Py{io) as a function of frequency inside and out- 
side the MNP. Blue circles (orange crosses) show a grid- 
ding of a/20 (a/ 10). The integration of the homogeneous 
GF over a cube of size a/20 (a/10) are given by blue-dark 
(orange- light) lines in (a-b). (a) z/a = in log scale, (b) 
z/a — 0.9 in log scale, (c) z/a = 1.2 in linear scale; here, 
the green solid line shows the analytic results. 

imply poor numerical convergence. Outside the MNP, we 
see that both FDTD grids give excellent agreement with 
each other and with the GF calculated using an analyti- 
cal scattering technique [10]. This confirms FDTD com- 
putes the correct regularized GF for a finite-size dipole. 
In Fig. 2 we examine the LDOS as a function of fre- 
quency corresponding to the scenario shown in Fig. 1 (b), 
for three different heights. Fig. 2(a): z/a = 0, Fig. 2(b): 
z/a = 0.9, and Fig. 2(c): z/a = 1.2 [indicated by the 
dashed lines in Fig. 1(e)]. We plot FDTD results in each 
as blue circles for a gridding of a/20, and orange crosses 
for a gridding of a/10. Additionally, we plot the regu- 
larized homogeneous GF calculated by integrating the 
homogeneous GF over a cubic volume of the exact same 
size as the FDTD grid, which we show as the solid lines 
in Fig. 2(a) and Fig. 2(b). We recognize that at the cen- 
ter of the MNP [Fig. 2(a)], the LDOS is dominated by 
the homogeneous contribution as there is no deviation 
between the homogeneous case and the total MNP case. 
In both cases the large peak at 3.23 eV corresponds to 
Re[e] = 0. We also see that the formal integration of the 
GF agrees perfectly (to within less than a percent) for 
both grids demonstrating again that the grid-dependent 
LDOS is actually the result of a finite sized dipole. As we 
approach the surface of the MNP, we see deviations from 
the homogeneous solutions [e.g., see region near 2.9 eV 
in Fig. 2(b)]; these deviations are due to the many non- 
dipolar surface plasmon modes that exist at the surface 
of the MNP [13] and are spectrally located in the same 
region as they are outside the MNP [Fig. 2 (c)]. In both 
Fig. 2 (a)-(b) we see that the peak of the LDOS is w lO'^ 
the value of the vacuum GF (for a/20), which indicates 
that the homogeneous GF would give spontaneous emis- 
sion enhancements (Purcell factors) which are orders of 
magnitude larger than any Purcell factors reported for 
semiconductor microcavities, though most of the emis- 
sion is into non-radiating modes. 
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Fig. 3: Local field calculations (real cavity) versus regu- 
larized GF calculations. In blue (right) we show the reg- 
ularized GF for the homogeneous case [circles, Fig. 1(a)] 
and for the total GF inside a MNP [crosses, Fig. 1(b)] 
using FDTD with a gridding of a/20. In orange (left) we 
show the local field GF for the homogeneous case [circles, 
Fig. 1(c)] and for the total GF inside a MNP [crosses. 
Fig. 1(d)] using FDTD with a gridding of a/20, where 
we have inserted a spherical cavity filled with vacuum 
where the cavity volume gives the same effective volume 
as an FDTD grid cell, (a/20)'^. The orange solid line cor- 
responds to exact calculations of the homogeneous local 
field GF using analytical techniques [Fig. 1(c)]. 

When we examine the LDOS outside the MNP, we 
observe that the difference between the two FDTD grids 
are much reduced compared with inside, and they agree 
very well with analytical techniques (green solid line). 
The peak Purcell factor here is < lO'' which is three or- 
ders of magnitude smaller than inside the MNP for the 
a/20 grid. We also recognize that at such short distances, 
the medium response is dominated by higher order plas- 
mon modes, and the dipolar localized surface plasmon 
mode at 2.77 eV is again a factor of 10 smaller. The clear 
agreement between FDTD with analytic techniques is 
encouraging for the calculation of spontaneous emission 
rates in arbitrary lossy material geometries. This is es- 
pecially true as we show FDTD is capable of calculating 
the GF both inside and outside a spherical MNP, which 
is one of the most difficult systems for FDTD given the 
underlying rectangular grid and the large index contrast. 

Next we examine local field effects to contrast with 
direct regularization of the GF, and to help describe the 
scenario of an embedded photon emitter with a different 
refractive index (e.g., a QD in a metal MNP). In Fig. 3 
we compare FDTD results for each of the four scenarios 
as shown in Figs. 1 (a)-(d) using a grid size of a/20. For 
the local-field cavity, we insert a small region of vacuum 
at the center of the MNP where the radius of the cav- 
ity yields the same effective volume as an FDTD grid 
cell, namely (a/20)'^. The homogeneous GF (circles) and 



the total GF (crosses) are shown for both the regularized 
case (blue-right) and the local field case (orange- left). We 
see that the local field case exhibits an increase, as well 
as a frequency shift when compared with the regularized 
GF, which is the result of the fundamentally different 
scattering geometry. To ensure these discrepancies are 
not caused by FDTD numerics, we calculate the homo- 
geneous local field case in Fig. 1(c) (solid orange line) 
using exact analytical scattering techniques [1, 10] and 
see very good agreement with the local field FDTD re- 
sults. Interestingly, in the local field case, the local field 
homogeneous case is still the dominant contribution just 
as in the regularized GF case and is the same for spher- 
ical or cuboid (not shown) cavities. For cavities of real 
£qd = 12 (typical for semiconductor QDs) the real cavity 
resonance retains the same height, but shifts to 2.24 eV. 

We summarize the regularization techniques as fol- 
lows. For the case where an emitter is located in a lossy 
material, but with no dielectric mismatch, then the regu- 
larized GF provides the proper method of regularization. 
However, for emitters with a dielectric mismatch such as 
QDs or atoms embedded in lossy host materials, the local 
field GF provides the correct manner of regularization. 
For the case when the QD itself is lossy, and is embed- 
ded inside a lossy material, then both regularization and 
local field effects must be considered. 

In conclusion, our numerical calculations presented 
here have shown that FDTD is an excellent tool for the 
calculation of spontaneous emission rates in and near 
lossy materials, and that the technique correctly takes 
into account regularization and local field effects. It is 
also one of the few tools available which allows the cal- 
culation of the GF in completely arbitrary geometries. 
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